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We present a calculation for the second moment of the local density of states in a model of a two- 
dimensional quantum dot array near the quantum Hall transition. The quantum dot array model 
is a realistic adaptation of the lattice model for the quantum Hall transition in the two-dimensional 
electron gas in an external magnetic field proposed by Ludwig, Fisher, Shankar and Grinstein. We 
make use of a Dirac fermion representation for the Green functions in the presence of fluctuations 
for the quantum dot energy levels. A saddle-point approximation yields non-perturbative results 
for the first and second moments of the local density of states, showing interesting fluctuation 
behaviour near the quantum Hall transition. To our knowledge we discuss here one of the first 
analytic characterizations of chaotic behaviour for a two-dimensional mesoscopic structure. The 
connection with possible experimental investigations of the local density of states in the quantum 
dot array structures (by means of NMR Knight-shift or single-electron-tunneling techniques) and 
our work is also established. 

PACS numbers: 71.55.Jv, 73.20.Dx, 73.40.Hm, 05.30.-d 



I. INTRODUCTION 



Recently, there has been a considerable surge of interest in the electronic properties of nanometric-scale metal 
and semiconductor structures |l],|| . Static as well as transport properties of these systems can be obtained from the 
knowledge of the statistics of the energy level distribution for quantum-mechanical structures in which weak disorder 
(or quantum chaos) plays a fundamental role. The presence of an external magnetic field enhances the influence of 
the quantum fluctuations on the properties of these mesoscopic systems. The remarkable advances in nanostructure 
fabrication procedures have also lead recently to a good deal of experimental and theoretical studies of the properties 
of Quantum Dots, Quantum Wires and, moreover, of Quantum Dot Arrays (QDA) || which can be obtained by a 
variety of techniques at semiconductor surfaces and interfaces. All these mesoscopic structures exhibit a wealth of 
new interesting quantum phenomena in an external magnetic field and at low temperatures. In particular, QDA (as 
"artificial crystals" ) offer the possibility of investigating situations not accessible in natural crystals for the magnetic 
fields that can be produced in normal laboratory conditions. 

On the other hand, a most celebrated quantum phenomenon, occurring at low temperatures in the presence of a 
magnetic field and assisted by the presence of weak disorder in a two-dimensional (2D) semiconductor heterostruc- 
ture, is the Integer Quantum Hall Effect (IQHE) Q. The plateaux in the Hall conductivity arc believed to be the 
result of current-carrying edge states and localized bulk states ||. The transitions between plateaux are due to 
localization-delocalization transitions induced by the combined effects of the magnetic field and disorder close to zero 
absolute temperature ||. The behaviour near a single IQHE transition (QHT) can be described in a number of ways, 
theoretically, in order to reproduce the jump in the conductivity and other singularities near the QHT. A model for 
the IQHE that is particularly suited to the present work has been proposed and studied by Ludwig et al. [j7|. In this 
model, a square lattice hopping model is defined, with a fixed half-flux-quantum per plaquette defining the lattice 
spacing a, with both nearest-neighbour and next-nearest-neighbour hopping processes considered in the presence of 
a staggered chemical potential and suitable forms of the one-particle disordered potential. In the continuum limit, 
a — > formally, and near the QHT this model yields a valid description of the behaviour near a single IQHE step 

§§■ 

It would seem natural, therefore, to extend the model of Ludwig et al. to the description of the QHT in a QDA 
system. This has indeed been done ||, by taking into account the fluctuations of the quantum levels pertinent to 
a single quantum dot. The existence of a supporting lattice in the QDA also makes the model of Ludwig et al., as 
extended and studied by Ziegler physically more appealing. In this model a QHT is indeed found and described 
in the proximity of the transition point. It is found that the otherwise sharp QHT step becomes rounded by the 
presence of level-fluctuations in the single-dots. 

In this paper we present an analytical characterization of the fluctuations of the local density of states (DOS) 
for a QDA near the QHT, using the model and calculations presented in M. Our calculation is interesting in its 
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own right for a good number of reasons. First and foremost, the averaged DOS and the moments of its distribution 
due to level statistics can be studied by means of Knight shift NMR measurements, now becoming accessible for 
semiconductor heterostructures [1C|. Another candidate technique for measuring local DOS-related properties is the 
single-electron tunneling (SET) 11], which provides an imaging of the local DOS in (bulk) semiconductors. A more 
detailed discussion will be presented in the Conclusions, Section IV. Also, from the more theoretical point of view, 
it is of interest to characterise the fluctuations near the QHT of an observable quantity like the DOS, considering 
that in the absence of fluctuations in the levels of a dot the DOS itself vanishes at the centre of the band, E = 0. 
When fluctuations are present, a semi-circle law is obtained for the averaged local DOS (p(r)) which does not vanish 
anymore at E = 0. The second moment of the global DOS vanishes identically, but not that of the local DOS 

M 2 (m) = (p(r,m) 2 )-( (0 (r,m)) 2 , (1) 

which we calculate in this work. Here, m plays the role of a chemical potential which controls the density of fcrmions 
in the system (for details see Section II). We find that Mi diverges near the point m = m c where (p(r,m)) vanishes, 
indicating that the level-fluctuations are particularly strong near this value of the energy. From the technical point of 
view, our work represents an attempt to characterise quantum mesoscopic fluctuations in a system of dimensionality 
D > 1 without applying the 2 + e expansion technique. Most approximate methods in fact deal with D = systems 
in the end (e.g. single dots, or aggregates of small metallic particles) and more rigorous and/or numerical work 
has been done chiefly for one-dimensional systems. An alternative approach to random electronic systems remains, 
naturally, the above-mentioned 2 + e expansion technique p^JT^ | , which however applies only to the scaling regime 
near the mobility edge. Therefore, our contribution to this field is instrumental for extending theoretical research to 
systems of dimensionality D = 2, as our calculations can indeed be implemented, with some additional computational 
effort, for the characterisation of the higher moments of the local DOS, as well as of those for the conductivity. As a 
further remark, we point out that the straight Quantum Hall system can be viewed as a mesoscopic-like system itself 
in terms of the recently-discovered jl4| universal conductance fluctuations between Quantum Hall plateaux near the 
QHT. This indicates, as confirmed by very recent numerical studies |ll|[l6|], that the characterization of the electronic 
states in the presence of level fluctuations for 2D electron-gas systems is of interest even before the complications of 
a periodic QDA structure are inserted in these devices. 

The article is organised as follows. In Section II we define the details of the model and justify its applicability to 
the QDA system in the close proximity of the QHT. This is a model for quasiparticles (quasielectrons) in a lattice of " 
artificial atoms" (the QDA); the quasiparticle levels are subject to some statistics with strong correlations of universal 
nature and characterized by a Wigner-Dyson distribution Jl7],[l8|]. The model is set up with a fixed flux i<& per 
plaquette, but although this may imply a fixed value of the magnetic field we stress that the presence of a (staggered) 
chemical potential /1 allows one to tune the model in the close vicinity of a QHT. The model includes a random 
potential, and we take the point of view that a random (Dirac) mass for the electrons is adequate in capturing the 
physics of the level fluctuations inside a single quantum dot. In Section III we briefly discuss the functional integral 
representation for the appropriate products of Green functions, following standard procedures [jl9|,^0). We set up a 
representation in terms of superfields (bosonic as well as fermionic) to take care of the averaging over level fluctuations, 
and carry out a large- N saddle-point approximation to evaluate the first and second moments of the local DOS near 
the QHT. This is a justified approximation when N w 10 2 represents the number of single quantum dot levels, but the 
case N = 1, representing the straight Quantum Hall system, can in principle also be described by this approximation. 
We carry out the analysis for the Gaussian fluctuations and present our results. These are discussed in the conclusions, 
Section IV, where the connection with some other recent papers on the subject is established, as is the connection 
between our results and their possible experimental verification by means of NMR Knight shift or SET measurements. 



II. THE MODEL 



A 2D array of quantum dots in a homogeneous perpendicular magnetic field can be modelled by a tight- 
binding tunneling (or hopping) Hamiltonian which, neglecting quasiparticle-interactions, takes the quadratic form 

H = J2 H"'" c r^ c r' + n - c -' w ^ n c r * creating a quasielectron at QDA site r in single-dot level a = 1, 2, . . . , N. The 

matrix elements H®'" are chosen as || 

H?*' = H^ a 'S ry + HW,5 a , a , + V r 5 r , r ,S a , a ,. (2) 

Here we take the point of view that the single-dot electronic states are statistically distributed with matrix-elements 
Hi- a defined by the Gaussian unitary ensemble (GUE) statistics: (Hr a ' a Hr ) = {9 / N)8 a .f35 a i ^1 , g being the 
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strength of the level fluctuations depending on the nature of the interactions and/or disorder and chaos inside a single 
dot. We take, at least initially, the point of view that weak tunneling processes take place between neighbouring dots 
and within the same single-dot energy levels. Strictly-speaking our formulation already contains tunneling between 
different neighbouring dots' energy levels via the random occupation of levels pertinent to each individual dot. Here, 
tunneling diagonal in the level indices will be assumed as realistic for the actual calculations and we take the tunneling 
rates to be non- vanishing only between nearest-neighbouring and next-nearest-neighbouring dots. These rates take 
into account the presence of a magnetic field that can be thought of as fixed at the value B = $o/2a 2 , in correspondence 
with the spacing a (typically H one has a = 100 to 500 nm) between the dots. In the Landau gauge, the tunneling 
Hamiltonian matrix elements read (setting r = (x, y) and indicating with e x and e y the lattice unit vectors): 

iJ«, = te™y/ a 6 r ,, r+esc + tS r , tr+ey ± it'e^ y ' a 8 r ,^ ±ey + h.c. (3) 

Finally, the potential V r represents an additional (e.g. electric) external field and can be regarded as a staggered 
chemical potential: V r = (— l) x+y /.i, which opens a gap 2fj, in the spectrum of the quasielectrons. For the "simple" 
IQHE situation, the model reduces to its N — 1 limit and corresponds to the model proposed for the QHT by Ludwig 
et al. 0. We assume for the QDA that the separation a between the dots is much greater than the dots size. 

Now, following a standard procedure @,|[, the model (||) is reduced in the continuum limit (formally a — > 0) and 
for long wavelengths to an equivalent 2D Dirac Hamiltonian [|]|| 

H al3 = ( ffl Vi + <7 2 V 2 )<Sa,/3 + <JsM a>3 (a,P = l,...,N). (4) 

V is here the 2D (lattice) gradient operator, <jj are Pauli matrices and M is a random mass matrix with mean 

m = n-t'\ = mSa^p + 6M?P with (8M^SM^') SM = (g/N)8 a ^5 a , ^5 r . r > . The tunneling rate t is scaled out 
in the Hamiltonian. Therefore, m is measured in units of t and the fluctuation parameter g is measured in units of t 2 . 
This model, for N = 1 and within a saddle-point approximation, describes the jump in the conductivity characterising 
the QHT and displays for the local DOS a characteristic semi-circle law ||p0|. The transition is driven either by 
changing the external field or, equivalently, as will be argued in this work, by changing the chemical potential fx. 

In principle there can be three different types of randomness: a random Dirac mass, a random energy and a random 
vector potential. We believe, however, that the random mass alone correctly accounts for the effects of the level- 
fluctuations in the single quantum dots. The random vector potential does not break the time reversal symmetry of 
the massless Dirac fermions. This is a strong restriction which implies that the case of the massless model has very 
special properties jl]]. Most remarkable is the fact that the average DOS is singular and obeys a power law: (p(E))^E a 
with — 1 < a < 1. This is an unphysical behaviour because the DOS should be finite at the Hall transition. A finite 
DOS exists for a random Dirac mass |],|2(J and for a random energy term ||l[]. (Ludwig et al. obtained for a 
random Dirac mass the result of a vanishing DOS at the Hall transition, using a perturbative renormalization group 
calculation. Consequently, they argued that only a combination of a random mass and a random vector potential can 
lead to a non-zero DOS. However, a non-perturbative approach gives a non-zero DOS already for a random Dirac 
mass proportional to expf— ir/g).) The random Dirac mass creates spontaneously a contribution proportional to 
the random energy (see Ref. || and Section III. A of this article). This implies that including the random energy will 
not lead to a qualitative change of the effect of the random Dirac mass. 

Results from previous calculations with this model will be recalled below; here it suffices to notice that the Green's 
function G = (H + ieoo) 1 has special properties because of the Lorentz invariance of the Dirac theory. In particular 
we will make use in the following of the relation (G"" rl .)* = — G^rr- Here, and in the following, oo is a convenient 
notation |22| for the 2x2 unit matrix. 

As a last remark about the model, we now clarify the role of the "fixed" magnetic flux per plaquette. The physical 
parameter of the model is not the magnetic flux, but the filling factor v = n$o/B (n is here the density of particles, B 
the external magnetic field) . n can be written as the total number of particles divided by the total area of the system. 
Since the particles are non-interacting fermions, there are at most N particles, if N is the number of lattice sites. 
The local DOS is symmetric around /i = 0. This implies that there are N/2 particles at /i = 0. On the other hand, 
the lattice can be divided into plaquettes of four lattice sites each. The flux through each plaquette is $0/2 = a 2 B, if 
a 2 is the area of the plaquette. Taking into account that there are N plaquettes (in the thermodynamic limit, when 
boundary effects can be ignored), we can write for the filling factor 

v = (N/2)<P /(Na 2 B) = (N/2)<S> /(N<I> /2) = 1. (5) 

Hence, the filling factor is 1, regardless of all other model parameters. 
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A. Second Moment of the Local DOS 



According to standard Green's function theory, the local DOS (LDOS) is obtained from 

1 



p<*<*(r,M) = -Im[G?? rr + G%? rr ] 



i 

2^ 



GQCt 



r~icxot 



V^ll.rr) l°"22, rr ) 



(6) 



For the tight-binding model without magnetic field this implies that one has to evaluate the two-particle Green's 
function with opposite signs of the frequency (advanced and retarded Green's functions) in order to get the second 
moment of the LDOS |l9| , |28| . This is the same problem as in the evaluation of the conductivity. In the model with 
half a flux quantum per plaquette, however, we can use the relation (GfFrr)* = — ^22 rr t° write 



p<*<*(r,M) = -(G^ rr + G^ rr ), 

7T 



and consequently 



p aa (r,M)p^(r,M) = - — (G^ r 

7T 



u 22.rr 



/3/3 



mi 



G 



0/3 



■)• 



(7) 



(8) 



In this way we set out to evaluate 

Mf = (p aa (r,M)p^(r,M)) 5M - ( P aa (r,M)) SM (p 0f3 (r,M)) SM 

= (jfi rr^ '11, rr) 'SM ~ (G\ l, rr )sM {G^i rr ) SM + (G" l, r r^22,rr) SM ~ (^ll,rr)sM {G^.rr) 



2]. (9) 



III. FUNCTIONAL INTEGRAL REPRESENTATION 



By means of a standard representation, we write, quite generally |20| 

XrjXrj exp(-Si) JJd$rd$r = -i(XrjX"j)s, 



G aa = — i 



and consequently 



J r 

with the supcrsymmetric action (sum convention for a) 

Si = -io- £ [($, (H + «r )*) + £ «5M r QQ '($^' ■ <r 3 $?)}, 



(10) 



(11) 



(12) 



where <7 e = sign(e) and the field = {^rj'Xrj)- The first component is Grassmann and the second complex. We 
notice that the normalization of the functional integral in ( |TT| ) is due to the combination of Grassmann and complex 
fields. Averaging with Gaussian distributed fluctuations in the "masses" yields exp(— S 2 ) = (exp(—Si))sM, with 



S2 = 



(13) 



Thus we have derived an effective field theory for $ which serves as a generating functional for the averaged Green's 
function. It is important to notice that not only SM creates the fermion-fermion interaction in (|T^ ) but that also 
other types of randomness can do this job. For instance, the interaction can also be created by a term which couples 
to a matrix field (fi = 1,...,4 includes the complex and Grassmann components): 



exp 



N 



exp 



with the supermatrix 
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Q e 
e -iP 



(15) 



The field $ appears only in a quadratic form on the right hand side. Therefore, it can be integrated out. This leads 
to 

(G^ tjj ) SM = J g njj exp(-NS 3 )V[Q, P, 9] = (Qujjfo, (16) 

and consequently 

(Grr,jjGrr,kk) SM = J (Sa,kjG2ijkS a ,fi + Gix^Q^kk) eM~NS 3 )V[Q, P, 9] 

= ( — Gl2,kjG21,jk5 a ,p+Gll,jjG22,kk)Q, (17) 

where we have defined 

[G q 1 -2tQt -2tQt y 1 

g -{ -2t9t G^+2 i rPr) rr - (18) 

Here G = (H Q + ie<7o) -1 , t = y^as, H = (H) = (ctiVi + (T2V2) + a 3 m and 

where we have introduced the graded (or Jl9| super-) trace and determinant. This implies for the local D9S 

±J2(p™(r,M)) SM = lEEWs = ^([rQrr]^. (20) 

a a j j 

For the correlation of the local D9S, it follows that 

^ 2 

= ^ ^(Gl2,kjG21,jk5 a ,f3 — GlX,jjG22,kk}Q + (Gil ,jj)c) (Gll,kk)( 



],k=l 



(21) 



A. Saddle Point Approximation 

The number of levels N appears in front of the action. Thus the effect of level fluctuations for N — * oo can be 
evaluated within a saddle point (SP) approximation p0|| . The SP equation reads 

S 1 

— [-TrQ 2 + logdet(G^ 1 - 2tQt)] = 0. (22) 

A second SP equation appears from the variation of P by replacing Q — > —iP. As an ansatz we take a uniform SP 
solution Qo = —iPo = —(l/2)[ir)<r3 + m s cro]. Then J22] ) leads to the conditions r\ = rjgl, m s = —mgl/(l+gl) with the 
integral I = J[{m + m s ) 2 + r] 2 + k 2 ]~ 1 d 2 k/2ir 2 . There is both a trivial and a nontrivial solution, 77 = and r\ 7^ 0. The 
latter is possible for m 2 < m 2 with m c — 2exp(— Tr/g). This means the level fluctuations shift the energy e — > e + r\ 
and the Dirac mass m — > m = m + m s , where rj(m) and m s (m) are solutions of the SP equations (p2|). The sign of 77 
is fixed by the condition that ij must be analytic in e; this leads to sign(ry) = sign(e). It is important to notice that 
r\ is proportional to the energy and not to the mass. That means, the random Dirac mass creates spontaneously a 
contribution which would be expected from a random energy term. This fact indicates that an additional random 
energy term may not qualitatively change the properties of the random Dirac mass. 
The average local D9S can now be directly calculated from ( |20| ) in SPA 

(p aa (r,M)) SM « = ^-Vm 2 -m 2 &(m 2 c - m 2 ), (23) 
ng 2irg 

where 9(x) is the Heaviside step function. We see that a semi-circle law is reproduced. In the following we will 
consider only the regime where the average local D9S is nonzero, i.e., m 2 < m 2 . 
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B. Gaussian Fluctuations 



In order to evaluate the second moment of the local DOS the Gaussian fluctuations around the SP must be 
calculated. Since Q, P and O are 2x2 matrices, the fluctuations can also be parametrised as 4-component vector 
fields: q% = SQu, q 2 = (SQi 2 + 6Q 2 i)/2, q 3 = -i(5Qi 2 - 5Q 2 \)/2, g 4 = 5Q 22 with analogous definitions for p x , ...,p 4 
and the Grassmann field ...,ip4 with t/j 2 = (O12 + ©2i)/2 and ip3 = — «(©12 — ©2i)/2. The action of the Gaussian 
fluctuations reads, in Fourier representation |2(J 

S « J ^2 0-k)^'{qk,^q-k,tJ.' + Pk.iiP-k,^ +2-ijj k ^ip^ k4l >)d 2 k 1 (24) 

with the fluctuation matrix . The stability matrix reads 



W - rl 2 Z T T T , (25) 



hi 


h2 


Il3 


hi \ 




1 22 


1 23 


Ph2 




^23 


^33 


pha 




— T* 


— T* 
J 13 


In I 



with p = \ij ' yf and /1 — m/2 + irj. In particular, for a vanishing wave vector we have 

1/ 5 - M * 2 /2tt|a*| 2 1/0-1/2* 

'«"' = ">= I S O^/A/. I 

1/S-1/27T l/ 5 -/i 2 /27r|/i| 2 . 

Thus the Gaussian fluctuations are always massive except for the critical points m = ±m c . This reflects the discrete 
symmetry of the Hamiltonian ([|) (for a more detailed discussion of the symmetry properties see Ref. |2(J ) . For k' ~ 
we obtain 

In = 1/g - m* 2 /2tt|m| 2 



/1 



2 



7 



13 



2tt|m| ImI 
2tt|m| N 



I u = 1/g - 1/2tt + o(fc' 2 ) 



'23 



27T 

^22/33 = 2/g - 1/tt, (27) 



that is, I14 is real. 



C. Expansion of Q 

Retaining only terms up to the second order, we have 

(Gl2,kjG2l,jk8 a ,P — Gll,jjG 22 ,kkjQ + (Gll,jj)Q{Gll,kk)Q ~ ^ ai ^{[GoT@TGo]rr,kj [GoT&TGo]rr,jk) 5Q, (28) 

with 

r _ ( {Ho - 2tQ t + iujao)- 1 \ , , 

y °~{ {H -2rQ r + iu;a )- 1 )- [ZJ) 

Now we have to perform the Grassmann integrations, ending up with 

(Qj'j",r'Qk'k",r")sQ = Ujij>i tV Uk>k>',v'{ll>v,r>'ipv',r")sQ = (1 /2)Uf jtt ,uUy k „ yl^'r" ,vv> > ( 30 ) 

with the unitary matrix 



G 



i o o cr 

J_ I o 1 -i o 

~ V2 \ 1 1 

,0 1, 



(31) 



The stability matrix ( p6| ) has one zero eigenvalue if fi is real, i.e. for m — ±m c . Consequently, there is a singularity 
in the Fourier components of I(fc) at k = 0. The leading behavior in the k ~ asymptotics of I -1 can be extracted as 



(e u en) ffg (Qxxe 22 )s Q \ _ i fir/ ir^ 
(e 22 eii) 4 Q (e 22 e 22 ) 5Q J 2 1 17? I44 1 / 1 



-1/4 



2tt + e 2ix 5 2tt - g 



cos{2x) + (g/2Tr\ f i\ 2 )cos(x)k 2 \ 2tt - g -2n + e^g } ' 

(32) 

where /i* = |/i|e 4:E . We consider weak fluctuations: higher order contributions in g in the numerator and in the 
coefficient of k 2 have been neglected. We make use of the following approximation 



{H - 2tQ t + iw<T )zJ w <5 r ,r' / (H 2 + fc 2 )^ 1 rf 2 fc/(2 7 r) 



'o _°J =(l/47r)^log(l + A 2 H- 2 ) 



/i* 

-/! 



(33) 



where A is a UV cut-off. For simplicity, we assume in the following A = 1. Then the second moment of the local DOS 
reads 



Mf=^" 2 E«^ / K u ,(k)d 2 k/(2nf 
l,V J 



(34) 



with #1,1 = I x l, Ki <2 = I 14 \ #2,1 = I41 1 and #2,2 = 1^ and 71 = //* log(l + \fi\ 2 ), 72 = -7^. Furthermore, 



M. 



2Re^{2ir - e 2lx g) - 2| A i| 4 (— 2tt + g) C [1 - cos(2a;) + (g cos(x)/2ir\Li\ 2 )k 2 }- 1 d 2 k/(2ir) 2 1 (35) 



(36) 



where C — [(1/4tt) log(l + \fi\ )] /2. The factor in front of the integral reads also 

7r~ 2 j47r[cos(4x) + 1] - 2 5 [cos(2x) + 1] j^ 4 . 
Integration with respect to k gives 

[1 - cos(2.t) + (g cos(x)/27r\fi\ 2 )k 2 ]^ 1 kdk/2'K ■ 



I Ml 



log{ r 



1 — cos(2a;) 



Then 



M. 



a(3 



-S a .aC 



2gcos(a;) 1 1 — cos(2:r) + (<7COs(x)/27r|/^| 2 ) 
4|/Lt| 6 [cos(4a;) + 1] r 1 - cos(2z) 



7T U ) J 



g cos(x) 



log{ T 



— cos(2a;) + (g cos(a;)/27r|^| 2 ) 



n\u\ 2 )}- 



(38) 



IV. DISCUSSION 



Making use of the evaluation described above, we can now display and discuss the dependence of M 2 on the effective 
chemical potential m = /1 — t' and on the amount of chaos or disorder in each single dot, g. 

According to the average DOS (|2^) there is a narrow band of bulk electronic states. Moreover, there are also edge 
states which carry the Hall current for m > m c . In the following we will discuss only the bulk states. At the band 
center the second moment of the LDOS is — 5 a p(m c /2g) 4 /it. The ratio \[M^I P at m = then is \/livn c j2g. 

Fig. 1 shows the dependence of M 2 (suitably scaled) on the parameters m and g. This figure, obtained for typical 
values of the parameters, displays an interesting dramatic increase as a function of m for the second moment M 2 of the 
local DOS as the edge of the energy band, m = m c , is approached. This (logarithmic) divergence of M 2 at the edge 
of the band indicates a power law singularity for the fluctuations near the QHT. This is in agreement with the results 
of the 2 + e expansion 25 2fj]. Also, notice that, as they should, the fluctuations decrease as randomness decreases, 
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except in the vicinity of m = m c . To our knowledge, these results, obtained for a realistic lattice of "artificial atoms", 
have not been discussed before (but see, however, below). The data in Fig. 1 show quite clearly the presence of 
important fluctuations in the local DOS as a consequence of single-dot level fluctuations, the fluctuations becoming 
perhaps critical at the edge of the band. 

The level fluctuations for an isolated quantum dot are described by a single random matrix 03 M"' 3 which remains 
from the Hamiltonian of Eq.(|j). The random matrix has, in contrast to the full Hamiltonian (|j), a continuous unitary 
symmetry. This generates in the N — > 00 limit a saddle point manifold which we have to integrate out. The remaining 
functional is a zero-dimensional nonlinear sigma model for the evaluation of the fluctuations of the DOS. The 
cross-over from the fully coupled quantum dots to isolated dots can be achieved by decreasing the tunneling rates t 
and t' . Since t was scaled out in the Hamiltonian, we must study for this purpose the rescaled mass m/t and the 
rescaled disorder parameter g/t 2 . Sending t ~ * 0, m c goes like 2t/ ^/exp(27rf 2 / g) — 1 — > y / 2g/-K. Therefore, the level 
fluctuations create a broad "band" which is surpressed by the tunneling term to a much narrower band. Moreover, 
the fluctuations of the local DOS are also reduced by the tunneling between the quantum dots. 

Notice that the strength of the fluctuations in the local DOS depends on the parameter m = fx — t' . This requires 
tuning the value of the chemical potential fi or of the next- near- neighbour hopping t' . Both could be achieved in 
experimentally realizable QDA by means of a varying gate voltage or a suitable design in the gate structure and in 
its voltages. 

It is important to discuss, at this point, the relevance of our calculation for perspective experimental observations 
on QDA. To our knowledge the local DOS can be accessed experimentally by two possibilities: the Knight shift of 
the NMR (nuclear magnetic resonance) lineshape due to the electron polarization |27j] and the SET (single electron 
tunneling) technique JTl| |. From these measurements one can gain access to the statistics of the local DOS mesoscopic 
fluctuations. 



A. NMR Knight shift spectroscopy 

We begin by discussing the NMR Knight shift technique, which has recently witnessed considerable revival both 
for mesoscopic semiconductor structures and for high-temperature superconductors. For zero-dimensional mesoscopic 
structures (single dots, small metallic particles or random aggregates of these) the Knight shift method has been 
discussed by Efetov and Prigodin (l^^^] and by Fal'ko and Efetov p{J]. For the one-dimensional disordered metal a 
related discussion for the connection between the local DOS and the Knight shift can be found in the work of Altshuler 
and Prigodin J3(| . Here we give a formulation for the Knight shift in a 2D QDA. We adapt the treatment for the Knight 
shift in a traditional solid (like, e.g. Li) [E7J to the situation in which a (2D) solid of "artificial atoms" is considered. 
Assume that a tight-binding description of quasielectron tunneling (or hopping) is appropriate (as advocated in the 
calculation of the previous Sections) and that a single-dot confining potential approximately of the square-well type 
can be used to describe the single "atomic" levels {E a , ip a (r)} of a single dot. Then a "band-structure" approach 
becomes a realistic calculation of the energy levels of the entire 2D QDA made up of these "atoms" and we can envisage 
the creation of "bands" of energy states from each "atomic" level E a . We label the 2D QDA energy levels, therefore, 
E a (k) — E a + J2 S e lk <5 and these will normally overlap due to the closeness of the levels E a of the single-dots. By 
building the many-electron wave-function in the normal way, the hyperfine interaction Hamiltonian can be averaged 
over the electronic degrees of freedom to give the contribution to the nuclear Hamiltonian arising from the electron 
polarization in the QDA that ultimately leads to the Knight shift formula 



Au(r) _ 8tt 
uq 3 



de X (e)p(r,e)(\u e (r)\ 2 )\ Ea{ii)=e , (39) 



at the nuclear site r. In the above formula, beside the local DOS which is the subject of the present investiga- 
tion, appear also the single-dot wavefunction u £ (r) (averaged over an equal-energy surface) and the electron spin 
susceptibility 

X(e) = ( /(e, -i) - f(e, „ ). ,x , (40) 



where the approximation holds at low temperatures (where the Fermi distribution function f(e, m s ) can be expanded). 
Since the low-temperature form of x(e) is delta-like, we essentially get that the Knight shift is proportional to the 
local DOS, evaluated at the Fermi level, however through a multiplicative constant that is in principle also randomly 
distributed, in a way that is beyond the scope of the present treatment. We conclude with the form 
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-^^-(K(r)\")\^ EFX (EF)p(E F ,r), (41) 
w 3 

for the expected Knight shift from a QDA. If we assume that the wavefunction at the Fermi level has much smaller 
fluctuations than the local DOS, then the moments of the local DOS can be extracted by looking at the moments of 
the Knight shift itself. 



B. Single-electron tunneling spectroscopy 



SET spectroscopy is well-known as a tool for the investigation of discrete single- and few-electron energy states in 
quantum dots fl3l| . Lerner and Raikh [ p2[ have proposed to apply the technique to study the weak-localization effects 
in disordered heterostructures, and the use of SET spectroscopy for imaging the local DOS in a disordered (bulk) 
semiconductor heterostructure has been recently reported [jll). Thus, SET spectroscopy appears to be a promising 
tool for an experimental investigation of local DOS fluctuations in a QDA. 

Lerner and Raikh fl32| show that, for the resonant-tunneling junction geometry, the measurement of the conductance 
Gt gives direct imaging of the local DOS, p: 

7re 2 5T 

Gt = ^mT p ' (42) 

with To, k and S geometrical parameters. A more complicated treatment ensues for a disordered material in the 
junction's barrier. The conductance presents mesoscopic fluctuations at low temperatures and one can show |32| that 
the second moments of local DOS and conductance are related by the formula 

M 2 (G t ) M 2 (p) 
(Gt) (p) 

However, Schmidt et al. pT| ] propose a double-barrier tunneling experimental geometry which leads through the 
measurement of the I — V characteristic to the relation 



6G(r) = 5{dI/dV) cx S(dp(r, E) / dE) , (44) 

for the local fluctuations of the conductance G of the disordered sample. Thus, SET spectroscopy can provide 
information, in principle, to local DOS statistical properties perhaps also in the case of a QDA setup. 

We mention, before concluding this subsection, the possibility of imaging the local DOS for surface electronic states 
also by means of scanning-tunneling microscopy (STM) j33|. This technique also relies on the measurement of the 
differential conductance, dl/dV, at local sample's sites; from this, the shape of the local DOS and its fluctuations' 
statistics can be reconstructed. This technique is perhaps more suitable for 2D systems like our QDA. 

Whether from the Knight shift in NMR or from the conductance in SET experiments, the measurement of the value 
of the local DOS in a real QDA should allow (via, e.g. repeated measurements) for a statistics of its fluctuations and 
the experimental determination of M%. This will be an indirect observation of the existence of level-fluctuations in 
single quantum dots. 



C. Beyond 2D: d = 2 + e, 1 and dimensions 



In dimensions higher than d = 2, there already exist calculations of the second moment of the local DOS for the 
situation in which randomness creates a mobility edge in the Anderson transition for a disordered metal |34f| . These 
calculations are based on 2 + e type expansions on the Q-matrix non-linear a model and lead to the result that 
Mi ~ (E — Ec)~^ 2 with /i2 = 2 + O(e). This shows the existence of the possibility of a divergence for the second 
moment of the local DOS near the mobility edge. Although the Anderson transition situation is not exactly the same 
QDA problem treated in the present paper, nevertheless these calculations show that a divergence of M 2 is to be 
expected. 

In lower dimensions, we recall the work of Altshuler and Prigodin, for d — 1 |pQ| , and of Efetov and Prigodin, for 
d = (2^] (a single dot in our picture), for the case of a disordered metal. In d = 1, exact diagrammatic techniques 
can be employed to work out the entire probability distribution for the local DOS of a weakly disordered metal. The 
resulting distribution depends strongly on the type (open or closed) of boundary conditions. Finally, in d = the 
non-linear-cr-model calculations involving supersymmetry of Ref. Jl9|,p8| can be employed. 







D. Conclusions 



We have shown how, by means of the technique of Dirac fermions in d — 2, one can characterize the distribution 
of the fluctuations of the local DOS in the proximity of an IQHE transition. We have restricted the calculation to 
characterize the second moment M%(p), but in principle also the higher moments can be evaluated (with increasing 
technical difficulties). 

We have made extensive use of the fact that in the Dirac fermions approach one is dealing with a unique saddle 
point solution, unlike in the case of the non-linear-cr-model approach. This is a consequence of the discrete symmetry 
of the model of Ludwig et al. |t]] for the case of a random mass, which we use for our calculations. This simple 
structure of the saddle point solution made the 2D calculation possible. 
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FIGURE CAPTIONS 

Figure 1. Second moment M 2 of the local density of states (LDOS) as function of effective chemical potential, m 
and strength of level fluctuations, g. 
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